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Abstract. 

The EGRET instrument on the Compton Gamma Ray Observatory (1991-2000) has 
positively detected high energy 7-ray emission from more than 67 active galaxies of 
the blazar class. The majority of the EGRET blazars are flat-spectrum radio quasars, 
which are characterized by inferred isotropic luminosities often as high as 3 x 10 49 
ergs s -1 . The remainder are BL Lac objects, some of which have been detected at 
TeV energies (> 250 GeV) by ground-based atmospheric Cherenkov telescopes. One 
of the remarkable characteristics observed in blazars is that the 7-ray luminosity often 
dominates the bolometric power in these sources. The detection of blazars by EGRET 
has undoubtedly been one of the highlights of the mission, and has forever impacted 
our understanding of the emission mcchanisms in these objects. In this article, we 
summarize the properties of EGRET blazars, and review the constraints that the 
EGRET observations place on the various models of 7-ray production in these sources. 



INTRODUCTION 

The earliest surveys of the sky above 100 MeV were carried out by the COS-B 
satellite, which detected 7-rays from 3C 273 [1], the first extragalactic 7-ray source. 
At the time EGRET was launched in 1991, 3C 273 was the only extragalactic source 
known to emit 7-rays. The subject of 7-ray blazars got its first dramatic boost 
after the EGRET detection of several high energy 7-ray blazars [2-4]. EGRET 
discovered two remarkable characteristics of these sources: (1) In the majority of 
the blazars, the 7-ray emission dominates the bolometric power in the spectral 
energy distributions of these objects, and (2) blazars are characterized by short- 
time scale variabilities, often on the order of days to months. 

Almost all the active galaxies detected by EGRET belong to the blazar-class of 
objects. Although the term "blazar" is not well-defined, it includes OVV quasars, 
BL Lac objects, and high polarization quasars (HPQs). These sources are char- 
acterized by emissions that include high radio and optical polarizations, strong 
variability at all wavelengths, and non-thermal, continuum spectra. Blazars de- 
tected by EGRET are radio-loud, flat-spectrum radio sources, with radio spectral 
indices a r > —0.6 for a flux density of S r oc z/* r . Several blazars have been ob- 



served to exhibit apparent superluminal motion (3C 279, 3C 273, 3C 454.3, PKS 
0528+134, for example [5]), as evidenced from VLBI radio observations. EGRET 
observations of blazars are summarized in section I. 

Three of the EGRET blazars (Mrk 421, Mrk 501, PKS 2155-304) have now been 
detected at TeV energies (> 250 GeV) by ground-based atmospheric Cherenkov 
telescopes (ACTs) [6] . These results have reinvigorated the 7-ray studies of blazars 
and have opened a new chapter in high energy astrophysics. As new ACTs with 
lower energy thresholds (< 50 GeV ) become operational [7], the high energy 7-ray 
studies of blazars will be conducted from the ground before the next generation 
satellite-based experiments like GLAST [8] or AGILE [9] come on line. The future 
of high energy 7-ray blazar studies looks very promising. Some of the EGRET 
blazars observed at TeV energies are described in section II. 

EGRET observations of blazars have helped constrain models of 7-ray emis- 
sion mechanisms in these objects. The combination of high luminosity, short-term 
variability, and superluminal motion observed in blazars all bolster arguments for 
relativistic beaming in these objects. This is in agreement with the strong asso- 
ciation of EGRET blazars with radio-loud flat-spectrum radio sources, with many 
of them showing superluminal motion in their jets. Most blazar models adopt a 
beaming paradigm. EGRET observations of blazars have enabled us to test several 
competing blazar emission models, most of which can be divided into two broad 
classes: "leptonic" and "hadronic." It is generally believed that the broadband 
spectral energy distribution of blazars is due to energy losses of highly energetic 
particles via both synchrotron and Compton processes. A summary of the current 
blazar models and their implications, as pertaining to EGRET sources, is given in 
section III. The possible contributions of blazars to the diffuse extragalactic 7-ray 
background are discussed in section IV. 

This article presents a snap shot of the current knowledge on EGRET-detected 7- 
ray blazars and selectively reviews some of the interesting highlights of the field. A 
comprehensive review is not possible, owing to the limited number of the allotted 
pages. Instead, the reader is referred to several recent articles which serve as 
excellent introductions to the exciting field of 7-ray blazars detected by EGRET 
[2,3,10,11]. 

I EGRET OBSERVATIONS 

This article reviews the blazars listed in the third EGRET (3EG) catalog [4]. 
These blazars were typically observed by EGRET for a period of about 2 weeks. 
Some observations were as short as 1 week, while others were as long as 3 to 5 weeks. 
EGRET's threshold sensitivity (> 100 MeV) for a single 2-week observation was 
~ 3 x 10 -7 photons cm -2 s" 1 . Details of the EGRET instrument, and data analysis 
techniques are given elsewhere [12,4]. 

Of the 271 sources detected by EGRET above 100 MeV, 67 are blazars detected 
with a high degree of confidence (a detection of at least a 4<r statistical significance 



at high galactic latitudes, and a 5a detection for |6| < 10° [4]). In addition, 27 other 
blazars are classified in the 3EG catalog as possible AGN [4] . This is because their 
identification is questionable, either because the candidate source has a low radio 
flux, or lies outside the 95% error contour of the EGRET source. This makes blazars 
the second-largest class of EGRET sources, following the "unidentified" sources 
that have no known counterparts at other wavelengths. EGRET has detected 
blazars with redshifts ranging from ~ 0.03 to 2.28. 

The majority of the blazars detected by EGRET are classified as flat-spectrum 
radio quasars (FSRQs). The 46 FSRQs with strong identifications, listed in the 
3EG catalog, are shown in Table 1. Of these, 0440-003, CTA 102, and 3C 454.3 
are not strictly "flat-spectrum" in the range between 1.4 and 5 GHz, or 2.7 and 8 
GHz. The BL Lac objects detected by EGRET are listed in Table 2. Not included 
in the table is Mrk 501, which is discussed in section II. The BL Lac objects are 
characterized by stronger polarization and weaker optical lines than FSRQs, which 
makes redshift determination more difficult. Other differences between FSRQs and 
BL Lac objects are discussed in the following sections. Not listed in either table 
are four flat-spectrum radio sources 0616-116, 0738+545, 1759-396, and 1908-201, 
also detected by EGRET. In addition, EGRET has detected 7-ray emission from 
Cen A, a radio galaxy, the only one seen to emit 7-rays [13]. 

A Time Variability of EGRET Blazars 

Blazars detected by EGRET have exhibited time variability in their 7-ray emis- 
sion on timescales of months, to sometimes well under a day. Long term flux 
variability for EGRET blazars have been discussed earlier by several authors [2,3]. 
Figure 1 shows the flux history of the blazar PKS 0528+134 over a period of six 
years [14]. Each data point corresponds to the time averaged flux from the source 
during a viewing period (VP). Upper limits, at the 95% confidence level, are shown 
for any detections below the 2a level. The figure demonstrates variability on the 
time scale of months. PKS 0528+134 is one of the most variable sources detected 
by EGRET. A dramatic flare was seen from the source in 1993 March when its flux 
was several times greater than that of the Crab [14]. 

The flux variability of the blazars can be quantified by assigning a "variability 
index," V = logQ, where Q = 1 — P x (x 2 , v) [15]. The quantity P x is the probability 
of observing x 2 or something larger from a x 2 distribution with v degrees of freedom. 
The x 2 values were obtained by fitting the data similar to that shown in Figure 1 
to a constant flux, using a least square fit method. The quantity V has been 
used in the past to assess the strength of variability [15,3], following the somewhat 
arbitrary classification: non-variable - V < 0.5; variable - V > 1; uncertain - 
0.5 < V < 1. Figure 1 shows a distribution of the variability indices of the FSRQs 
and BL Lac objects listed in tables 1 and 2. Although the sample of the blazars is 
small, the figure indicates that FSRQs are generally more variable than the BL Lac 
objects, in the EGRET energy range. Nearly 75% of all FSRQs may be classified 



TABLE 1. Flat-Spectrum Radio Quasars (FSRQs) detected by EGRET 



Name 


Other Name 


l 


b 


FluxxlO" 8 
ph cm~ 2 s 


Photon 
Spectral Index 




z 


L48 
ergs s -1 


0202+149 


4C+15.05 


147.95 


-44.32 


23.6 + 5.6 


2 


.23 ± 0. 


.28 





.405 





.0600 


0208-512 


PKS 


276.10 


-61.89 


85.5 + 4.5 


1 


.99 ± 0. 


,05 


1 


.003 


1 


.4500 


0336-019 


CTA026 


188.40 


-42.47 


118.8 + 22.0 


1 


.84 + 0. 


,25 





.852 


1 


.4000 


0414-189 




213.90 


-43.29 


49.5+16.1 


3 


.25 + 0. 


,68 


1 


.536 


2 


.1600 


0420-014 




194.88 


-33.12 


50.2+10.4 


2 


.44 + 0. 


,19 





.915 





.6900 


0430+285 




170.48 


-12.58 


22.0 + 2.8 


1 


.90 + 0. 


,10 




. 






0440-003 


NRAO190 


197.39 


-28.68 


79.0+10.1 


2 


.37 + 0. 


,18 





.844 





.9100 


0446+112 




187.86 


-20.62 


109.5 + 19.4 


2 


.27 + 0. 


,16 


1 


.207 


2 


.7900 


0454-463 




252.40 


-38.40 


7.7 + 2.1 


2 


.75 ± 0. 


,35 





.858 





.0900 


0459+060 




193.99 


-21.66 


12.1 + 3.1 


2 


.36 ± 0. 


,40 


1 


.106 





.2500 


0458-020 




201.35 


-25.47 


11.2 + 2.3 


2 


.45 ± 0. 


,27 


2 


.286 


1 


.1900 


0528+134 




191.50 


-11.09 


93.5 + 3.6 


2 


.46 + 0. 


,04 


2 


.060 


7. 


.8600 


0827+243 




199.91 


31.69 


24.9 + 3.9 


2 


.42 ± 0. 


,21 


2 


.046 


2 


.0600 


0836+710 




143.49 


34.79 


10.2 + 1.8 


2 


.62 ± 0. 


16 


2 


.172 





.9700 


0954+556 




159.55 


47.33 


9.1 ± 1.6 


2 


.12 ± 0. 


,18 





.901 





.1200 


1156+295 


4C+29.45 


201.53 


78.63 


50.9 + 11.9 


1 


.98 ± 0. 


.22 





.729 





.4300 


1222+216 




254.91 


81.53 


13.9 + 1.8 


2 


.28 ± 0. 


.13 





.435 





.0400 


1226+023 


3C273 


289.84 


64.47 


15.4+1.8 


2 


.58 ± 0. 


.09 





.158 





.0100 


1229-021 




292.58 


59.66 


12.7 + 2.9 


2 


.85 + 0. 


.30 


1 


.045 





.2400 


1243-072 




300.96 


55.99 


9.8 + 2.1 


2. 


.73 ± 0. 


.17 


1 


.286 





.2900 


1253-055 


3C279 


304.98 


57.03 


179.7 + 6.7 


1 


.96 ± 0. 


.04 





.538 





.7800 


1331+170 




346.29 


76.68 


9.4 + 2.7 


2. 


.41 ± 0. 


.47 


2 


.084 





.8100 


1406-076 




334.23 


50.30 


97.6 + 9.1 


2 


.29 ± 0. 


.11 


1 


.494 


4 


.0000 


1424-418 




321.66 


16.98 


29.5 + 5.3 


2 


.13 ± 0. 


.21 


1 


.522 


1 


.2600 


1510-089 




351.49 


40.37 


18.0 + 3.8 


2 


.47 + 0. 


.21 





.361 





.0300 


1606+106 


4C+10.45 


23.51 


41.05 


34.9 + 5.6 


2 


.63 ± 0. 


.24 


1 


.226 





.9200 


1611+343 




55.44 


46.29 


26.5 + 4.0 


2 


.42 ± 0. 


.15 


1 


.401 





.9400 


1622-297 




348.67 


13.38 


258.9 + 15.3 


2. 


.07 + 0. 


.07 





.815 


2. 


.7700 


1622-253 




352.28 


16.37 


42.6 + 6.6 


2 


.21 ± 0. 


13 





.786 





.4200 


1633+382 


4C+38.41 


61.21 


42.26 


107.5 + 9.6 


2 


.15 ± 0. 


.09 


1 


.814 


6 


.8000 


1725+044 




27.27 


20.62 


17.9 + 4.1 


2 


.67 + 0. 


.26 





.296 





.0200 


1730-130 


NRAO530 


12.00 


10.57 


36.1 + 3.4 


2. 


.23 ± 0. 


.10 





.902 





.4800 


1739+522 




79.37 


32.05 


18.2 + 3.5 


2 


.42 ± 0. 


23 


1 


.375 





.6200 


1741-038 




22.19 


13.42 


21.9 + 5.3 


2 


.42 ± 0. 


.42 


1 


.054 





.4100 


1830-210 




11.92 


-5.50 


26.6 + 3.7 


2 


.59 + 0. 


.13 


1 


.000 





.4500 


1933-400 




358.65 


-25.23 


21.9 + 4.9 


2 


.86 + 0. 


.40 





.966 





.3400 


1936-155 




23.95 


-17.12 


55.0 + 18.6 


3 


.45 + 1. 


.27 


1 


.657 


2 


.8400 


2022-077 




36.72 


-24.40 


74.5 ± 13.4 


2. 


.38 + 0. 


.17 


1 


.388 


2 


.5900 


2052-474 




352.56 


-40.20 


23.6 + 6.0 


2. 


.04 + 0. 


.35 


1 


.489 





.9600 


2209+236 




81.83 


-25.65 


14.6 + 4.2 


2. 


.48 ± 0. 


.50 










2230+114 


CTA102 


77.45 


-38.50 


19.2 + 2.8 


2 


.45 ± 0. 


.14 


1 


.037 





.3500 


2251+158 


3C454.3 


86.05 


-38.30 


53.7 + 4.0 


2 


.21 + 0. 


.06 





.859 


0. 


.6500 


2320-035 




76.82 


-58.07 


38.2 + 10.1 








1 


.411 


1. 


.3800 


2351+456 




113.39 


-15.82 


14.3 + 3.7 


2. 


.38 + 0. 


.38 


1 


.992 


1. 


.1200 


2356+196 




107.01 


-40.58 


16.0 + 4.7 


2. 


.09 + 0. 


.35 


1 


.066 





.3100 



TABLE 2. BL Lac Objects detected by EGRET 



Name 


Other Name 


l 


b 


FluxxlO" 8 


Photon 




z 


Las 










ph cm~ 2 s _1 


Spectral Index 






ergs s -1 


0219+428 


3C66A 


140.22 


-16.89 


18.7 ± 2.9 


2 01 ± 14 





.444 


0500 


02354-164 


OD+160 


156.46 


-39.28 


65 1 ± 8 8 


1 85 + 012 





.940 


9500 


0454-234 




223.96 


-34.98 


14.7 ± 4.2 


3.14 ± 0.47 


1 


.009 


0.2500 


0537-441 




250 08 


-30.86 


25.3 ± 3.1 


2.41 ± 0.12 


o 


.894 


3300 

KJ ■ KJ KJ KJ KJ 


0716+714 




143.98 


28.00 


17.8 ± 2.0 


2.19 ± 0.11 


o 


.300 


0200 

V'.V/-— V/V/ 


0735+178 

1 KJKj | i- 1 l_J 




202.16 


17.88 


16.4 ± 3.3 


2 60 ± 28 


o 


.424 


0.0400 


0829+046 




219.60 


23.82 


16.8 ± 5.1 


2 47 ± 40 


o 


.180 


0.0100 


0851+202 


OJ+287 


207.19 


35.43 


10.6 ±3.0 


2.03 ±0.35 





.306 


0.0100 


0954+658 




145.78 


43.11 


15.4 + 3.0 


2.08 + 0.24 





.368 


0.0300 


1101+384 


Mrk421 


179.97 


65.04 


13.9 + 1.8 


1.57 + 0.15 





.031 


0.0002 


1219+285 


ON+231 


197.27 


83.52 


11.5 + 1.8 


1.73 + 0.18 





.102 


0.0015 


1334-127 




320.07 


46.95 


11.8 + 3.4 


2.62 + 0.42 





.539 


0.0500 


1604+159 


4C+15.54 


29.18 


43.84 


42.0 ± 12.3 


2.06 + 0.41 





.357 


0.0800 


2032+107 




56.12 


-17.18 


13.3 ±3.1 


2.83 + 0.26 





.601 


0.0700 


2155-304 




17.45 


-52.23 


30.4 + 7.7 


2.35 ±0.26 





.116 


0.0100 


2200+420 


BLLac 


92.56 


-10.39 


39.9 ± 11.6 


2.60 ±0.28 





.069 


0.0024 



Note that the table lists sources only from the 3EG catalog, and does not include Mrk 501, 
which was not detected by EGRET until later. The EGRET detection of this source is 
described in section II. 



as variable. In comparison, most of the unidentified EGRET sources are found to 
be non-variable, although recent studies have found short time scale variability in 
a few unidentified EGRET sources [16]. 

Although the study of short time scale variability is limited for most blazars 
by the small number of photons detected by EGRET during a single observation 
period, some EGRET blazars have been observed to flare dramatically on the time 
scale of days. For example, PKS 1622-297 exhibited a major flare during the 
observation in 1995. A flux increase by a factor of at least 3.6 was observed in a 
period of less than 7 hours [17]. The peak flux observed was (17 ± 3) x 1CT 6 cm -2 
s -1 (E > 100 MeV), which corresponded to an isotropic luminosity of 2.9 x 10 49 ergs 
s -1 [17]. The characteristic time scale of variability, or the "doubling time," can be 
somewhat arbitrarily defined as the shortest timescale during which the source flux 
is seen to change by a factors of 2. In the case of PKS 1622-297, the doubling time 
was less than 3.8 hr. A dramatic flare was also observed in the blazar 3C 279, in 
which the doubling time was ~ 8 hours [18]. Other blazars that have shown strong 
variability on timescales of 1 to 3 days in the past are PKS 0528+134, 3C 454.3, 
PKS 1633±382, 1406-076, and CTA26. Studies of fast variations of 7-ray emission 
in blazars using structure function analysis techniques have also been carried out 
for a few selected sources [19,20]. 

The study of short time scale variability in blazars is important as they help to 
constrain the size of emission regions in these objects. A flux variation by a factor 
of two on an observed time scale St b s limits the size r of the emitting region to 
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EGRET Observations of PKS 0528+134 
E > 100 MeV 
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1. (Left) Flux history of the blazar PKS 0528+134. (Right) Distributions of variability 
FSRQs (dashed line) and BL Lac objects (solid line). 



roughly r < cSt bs/(l + z) using simple light travel time arguments. This implies 
that the blazar emission region is very compact. 



B Spectra of EGRET blazars 

The spectra of blazars detected by EGRET cover the energy range from 30 
MeV to 10 GeV and are well-described by a simple power-law model of the form 
F(E) = k(E/E Q )~ a photons cm" 2 s" 1 MeV" 1 , where the photon spectral index, a, 
and the coefficient, k, are the free parameters. The energy normalization factor, 
Eq, is chosen so that the statistical errors in the power law index and the overall 
normalization are uncorrelated. Details of the EGRET spectral analysis techniques 
for blazars may be found elsewhere [2]. Figure 2 shows the spectra of two blazars 
detected by EGRET [4]. The spectral indices of all the blazars detected by EGRET 
are included in tables 1 and 2. The average photon spectral index in the EGRET 
energy band, assuming a simple power law fit to the spectrum, is ~ 2.2. There is 
no evidence of a spectral cutoff for energies below 10 GeV. Figure 3 shows a plot 
of spectral indices versus redshift of all blazars detected by EGRET. 



C Luminosity of EGRET blazars 

Tables 1 and 2 list the 7-ray luminosity of the blazars, calculated assuming 
isotropic emission. The 7-ray luminosity was estimated using the method described 
in [3]. Figure 3 shows a plot of the luminosity as a function of redshift for the 
FSRQs and BL Lac objects detected by EGRET. The solid line corresponds to 
the typical detection threshold for EGRET as a function of redshift, for relatively 



good conditions. As seen from the figure, the BL Lac objects (indicated with dark 
diamonds) are predominantly closer and have lower luminosities than the FSRQs. 

II SEDS OF EGRET BLAZARS 

A FSRQs and Radio-selected BL Lac Objects 

In order to decipher the underlying physics of the emission mechanisms in blazars, 
it is essential to study the broad band spectral energy distributions (SEDs) of these 
sources. Since EGRET first started detecting blazars, several multiwavelength 
campaigns have been organized to simultaneously observe these sources at different 
wavelengths (see [21] for a recent review). Although these campaigns have proven 
difficult to organize, and the SEDs of blazars are often non-simultaneous or under- 
sampled, they have helped constrain blazars models. 

The SED of the blazar 3C 279 from radio to 7-ray energies over 15 decades in 
energy for two different epochs [18] is shown in Figure 4. 3C 279 is one of the 
most intensely monitored of all EGRET blazars. The figure plots the quantity vF v 
as a function of frequency, thus showing the relative amounts of power detected 

3EG_J1255-0549 3EG_J1222 + 2841 
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FIGURE 2. Photon spectra of the FSRQ 3C 279 (left) and the BL Lac object 1219+285 (right) 
in the energy range 30 MeV to 10 GeV. The solid lines are the best fit to a power law. 2a upper 
limits are shown as downward arrows. 



in each energy band. The overall SEDs of both FSRQs and BL Lac objects show 
two broad components: the synchrotron component and the inverse Compton (IC) 
component. In the case of FSRQs, as shown in Figure 4, the first peak is in the 
optical-IR band, while the second peak is in the MeV-GeV energy range. The figure 
demonstrates the dominance of the 7-ray luminosity over that at other wavebands. 
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FIGURE 3. (Top) Spectral index versus redshift for EGRET blazars. (Bottom) Luminosity vs 
redshift for blazars detected by EGRET. The BL Lac objects are indicated with filled symbols. 
The typical detection threshold for EGRET is shown as a solid curve. 




FIGURE 4. (Left) SED of 3C 279 from radio to 7-ray energies for diffcrcnt cpochs. The figure 
demonstrates the considerable variability that is often seen in the broad band spectra of blazars. 
(Right) Broadband SED of MrK 501 during several diffcrcnt cpochs. The filled and open squares 
correspond to Whipple flux on April 7 and 16, 1997, respectively. The circles correspond to data 
taken in 1996 March. The diamonds correspond to non-simultaneous data from the NED data 
base. (See text for references.) 

3C 279 also shows considerable spectral variability, particularly in the 7-ray band, 
between different epochs. This has been found to be a characteristic feature of the 
FSRQs observed by EGRET. 

SEDs of the FSRQs and BL Lac objects detected by EGRET have been studied 
by many researchers, both for individual sources as well as for source classes. Some 
recent work may be found in [22-25]. 

B Blazars detected at TeV Energies 

Of the blazars detected by EGRET, only three have been seen at TeV energies 
(> 250 GeV), all three being X-ray selected BL Lac (XBL) objects. These are 
Mrk 421, Mrk 501, and PKS 2155-304 [26]. Figure 4 shows the broad band SED 
of Mrk 501 measured during the outburst in 1997 May [27]. Mrk 501 was detected 
first at TeV energies before it was recognized in the EGRET data. The strongest 
EGRET detection of the source occurred in 1996 May, when Mrk 501 was detected 
at 5.2 o at energies above 500 MeV [28]. Six high energy photons (> 1 GeV) were 
observed from the direction of Mrk 501 in a 1-day interval, when the probability 
for this happening by chance was < 10~ 6 [28]. A power-law fit to the EGRET data 
suggested a spectral index of 1.2 ± 1.0 for Mrk 501, and is the hårdest known blazar 
spectrum at GeV energies [28]. 

The SED of Mrk 501 shown in Figure 4 is similar to that of other X-ray selected 
BL Lac objects. Unlike in FSRQs, the 7-ray peak in the SED of this high-frequency- 
peaked BL Lac (HBL) no longer dominates the spectrum, although it plays a 



very important part. The synchrotron peak of Mrk 501 is located in the X-ray 
band, while the Compton peak is in the GeV-TeV regime. In faet, a correlation 
between the peak frequencies of the synchrotron and IC components of the SED 
and the total energy density of the emitting region has been noted previously 
[24] . There appears to be evidence for a well-defmed sequence in the properties of 
different blazar classes, namely, HBL, LBL (low frequency-peaked BL Lacs), HPQ 
(high-polarization quasars), and LPQ (low-polarization quasars). The trends in 
the observed properties are a decrease in frequencies of the synchrotron and inverse 
Compton peaks, and an increase in the power-ratio of the high and low energy 
spectral components along the sequence HBL — > LBL — > HPQ — > LPQ. 

At a redshift of 0.031, Mrk 421 is the closest XBL seen by EGRET, and the first 
one to be detected by ground-based ACTs [6]. Although it is a rather weak source 
at GeV energies, it has been persistently detected by EGRET [29]. Above 30 MeV 
the source has a very hard photon spectral index of 1.57 ± 0.15, and has a high 
energy peak in its SED at TeV energies. 

III GAMMA RAY EMISSION MODELS 

A Leptonic models 

It is generally believed that blazars are powered by accretion of matter on to 
a supermassive black hole, and that 7-ray emission originates in strongly beamed 
jets. The absence of intrinsic 7 — 7 pair absorption in the observed blazar spectra 
strongly points to the faet that 7-ray emission is beamed radiation from the jet. A 
review of the constraints placed on AGN models from the 7-ray observations may 
be found in [30]. 

Leptonic jet models explain the radio to UV continuum from blazars as syn- 
chrotron radiation from high energy electrons in a relativistically outflowing jet 
which has been ejected from an accreting supermassive black hole [31]. The emis- 
sion in the MeV-GeV range is believed to be due to the inverse Compton scattering 
of low-energy photons by the same relativistic electrons in the jet. However, two 
main issues remain questionable: the source of the soft photons that are inverse 
Compton scattered, and the structure of the inner jet, which cannot be imaged di- 
rectly. The soft photons can originate as synchrotron emission from within the jet 
(the synchrotron-self- Compton or SSC process), or from a nearby accretion disk, or 
they can be disk radiation reprocessed in broad-emission-line clouds (the external 
radiation Compton or the ERC process). A review on leptonic jet models and the 
relevant references may be found in [32]. 

Recently, Mukherjee et al. (1999) have applied a leptonic jet model to the SED of 
PKS 0528+134 observed by EGRET from 1991 to 1997 [14]. A combination of SSC 
and ERC models was used. A two-component model, in which the target photons 
are produced externally to the 7-ray emitting region, but also including an SSC 
component, was required to suitably reproduce the spectral energy distributions of 




the source. The analysis indicated that there is a trend in the observed properties 
of PKS 0528+134, as the source goes from a 7-ray low state to a flaring state. The 
model fits presented in [14] indicated that during the higher 7-ray states, the bulk 
Lorentz factor of the jet increased and the ERC component dominated the high- 
energy emission. The energies of the electrons giving rise to the synchrotron peak 
decreased, and the power-ratio of the 7-ray and low energy spectral components 
increased, as the source went from a low to a high 7-ray state. 

In Figure 5 the fit results to two extreme states of PKS 0528+134 from [14] are 
compared. During the 1993 March observations, PKS 0528+134 was in its highest 
state ever, while the 1997 February data correspond to one of the lowest states 
of the source during the EGRET observations. The source exhibits considerable 
spectral variability between different epochs, a characteristic typical of many FS- 
RQs observed by EGRET. The model calculations of the SED of PKS 0528+134 
yielded a shift of the synchrotron peak frequency towards lower frequencies in a 7- 
ray bright state of the source, as has been predicted analytically [33]. The spectral 
variability in PKS 0528+134 appears to arise from the different Doppler boosting 
patterns of the SSC and the ERC radiations. The relative contributions of the 
SSC and ERC cooling mechanisms seem to be related to the optical to 7-ray flux 
ratio from the source. The SSC mechanism plays a larger role if the source is in a 
low flux state. The ERC mechanism is the dominant cooling mechanism when the 



source is in a high 7-ray state. 

Bottcher [33] has demonstrated that the predicted peak shift is characteristic 
of a leptonic jet model in which the 7-ray emission is dominated by the ERC 
process during 7-ray flares, but does not depend on the physical details of such a 
model. This is particularly important, as the SEDs of most blazars are still quite 
poorly constrained observationally, and can in many cases be fitted equally well 
with several fundamentally different models. 

A similar analysis was recently performed for the multi-epoch SEDs of 3C 279 
[34] , where the same general trend in the SEDs were observed as in the case of PKS 
0528+134. 

B Hadronic models 

In contrast to leptonic jet models where the jet consists of a hot e ± plasma, in 
hadronic models the observed 7-ray emission is initiated by hot protons interacting 
with ambient gas or low-frequency radiation. In the proton induced cascade (PIC) 
model [35,36] the relativistic protons interact electromagnetically and hadronically, 
i.e. through photo-pion production, with low frequency synchrotron radiation. 
Such a model would require extremely high luminosity in high energy protons in 
order to explain the observed TeV 7-ray fluxes in blazars. In contrast, in the 
synchrotron proton blazar models [37,38], the high energy emission in the blazars 
detected at TeV energies is attributed to synchrotron radiation of extremely high 
energy protons. In such models, synchrotron radiation becomes the dominant en- 
ergy loss mechanism of the high energy protons. Synchrotron proton blazar models 
have recently been applied to the case of the EGRET blazars Mrk 501 and Mrk 
421 [37,39]. 

IV BLAZARS AND THE EXTRAGALACTIC 
GAMMA-RAY BACKGROUND 

The precise origin of the extragalactic diffuse emission is not well-known and pos- 
sibly includes contributions from diffuse origins as well as unresolved point sources. 
Recently, the evolution and luminosity function of the EGRET blazars was used 
to estimate the contribution to the diffuse extragalactic background [40]. It was 
found that the evolution was consistent with a pure luminosity evolution. The 
redshift distribution of EGRET blazars was used to characterize the low end of 
the luminosity function better. The high end of the luminosity function was fixed 
by a non-parametric estimate [40]. Using the lower limit of the de-evolved lumi- 
nosity function, the 7-ray loud AGN contribution to the extragalactic 7-ray flux 
was estimated to be 4.0io;g x 10~ 6 photons cm -2 s -1 sr" 1 . The sky-averaged flux 
contribution of identified EGRET blazars is ~ 1 x 10 -6 photons cm" 2 s -1 sr _1 . 
Contribution to the diffuse background by unresolved blazars, therefore, was found 



to be ~ 3.0tg;g x 10 6 photons cm 2 s 1 sr 1 . The estimated extragalactic diffuse 
flux for E > 100 MeV is ~ 1.36 x 10~ 5 photons cm' 2 s' 1 sr" 1 [41]. These calcula- 
tions indicate that blazars cannot account for all of the diffuse extragalactic 7-ray 
background at the energies considered. Under this scenario, only ~ 25% of the 
diffuse extragalactic emission measured by SAS-2 and EGRET can be attributed 
to unresolved 7-ray blazars. These Åndings are consistent with other studies [42], 
and contrary to estimates which assume a linear correlation between the measured 
radio and 7-ray fluxes and attribute a much larger contribution by blazars to the 
7-ray background (e.g. [43]). These results lead to the exciting conclusion that 
other sources of diffuse extragalactic 7-ray emission must exist. 

V SUMMARY AND FUTURE PROSPECTS 

During its lifetime EGRET detected more than 67 blazars at energies above 100 
MeV. This has been one of the most remarkable contributions from the mission. 
Only a handful of the EGRET sources have been detected with ground-based ACTs 
at energies above 250 GeV, although most of them have been observed [44]. This 
indicates that spectral changes must exist in the unexplored range between 50 
and 250 GeV, possibly due to intrinsic absorption at the source or due to 7-7 
pair production off the extragalactic background light (EBL) [45]. In the future, it 
will be important to complement blazar studies by space-based experiments such as 
AGILE or GLAST, with ground-based instruments sensitive to the energy range 50 
to 250 GeV. Ground-based instruments like STACEE and CELESTE with energy 
thresholds lower than 250 GeV are currently operational and well-suited to the 
study of several northern hemisphere AGN. Other experiments like VERITAS, 
HESS, MAGIC and CANGAROO III will soon come online. A review of these 
experiments, and a few others, and the relevant references may be found in [7,26]. 
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